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To dear memory of Jiirgen Moser 



Abstract. A few remarks on integrable dynamical systems inspired by dis- 
cussions with Jiirgen Moser and by his work. 

1. Introduction 

Nine years have passed since the last time I spoke to Jiirgen Moser but it is 
stih painful for me to think that this will never happen again. That time it was in 
his house in Schwerzenbach near Zurich, which I was passing by returning from a 
conference in July 1999. It was after the metastases had been discovered in his body 
and he already knew his fate. He started, seeing a question in my eyes, by saying 
that he most probably had only up to Christmas time to live. Then he looked at 
the 1958 German record of St Matthew Passion I had brought to him and changed 
the subject by telling me about his time at Princeton, of his friendship with John 
Nash and many other interesting details of his life around 1958. It was the usual 
Moser I had been lucky to know for about 10 years. I am telling this because it 
shows how strong this man was, which in combination with his dignity and wisdom 
makes him unique to everybody who knew him well enough. 

I do not ever remember Moser arguing. It does not mean that he did not make 
clear his opinion, he actually always did this immediately in a very straightforward 
fashion but in such a way that would not offend even the most sensitive person. 
However he had never repeated his point afterwards in contrast, say, to me. On 
every occasion that I remember this happening, I was wrong and he was right. 
Once this was about the Euler rigid body system. I said that it seems to have 
been explicitly integrated in elliptic functions only in the second half of the XlX-th 
century (I think I mentioned Hermite in this relation). Moser did not think that it 
was that late but when I insisted he reacted in his usual way: "Really ?" The next 
day I found on my desk a library volume of Jacobi's Gesammelte Werke opened at 
the first page of Jacobi's letter (more than 50 pages long !) from March 1850 to the 
Academic des Sciences de Paris with full integration in elliptic functions not only 
of the Euler system but of the whole frame motion. There were other cases I had 
to understand my foolishness myself but Moser always left it with me to find out 
the truth and never returned to this point. I was not an exception of course: many 
of the visitors of FIM at ETH, Zurich during the time of his Directorate probably 
had similar experiences. 

My happiest memories of Moser are, when after his busy days, he invites me to 
his office to talk maths. He is sitting in his chair with glasses on and from time 
to time making notes. Taking notes was one of his habits, which the participants 
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of the Analysis seminar at ETH remember very well: he did this at every talk he 
attended. To talk to Moser was a real joy but sadly strict timekeeping was another 
habit of his, so to prolong these happy hours I sometimes accompanied him on his 
way to Stadelhofen railway station. 

This note, written for the special issue of "Regular and Chaotic Dynamics" 
dedicated to 80th anniversary of Moser, consists of a few remarks on integrable 
systems, which were inspired by discussions with Jiirgen Moser and by his work. I 
have presented them in several talks, but have never published them before. 

2. Integrable modulations: could the adiabatic invariants be exact ? 

One of the problems I discussed with Jiirgen Moser in the early 1990s, when we 
studied the integrable discretisations of the classical systems [M] , was the following. 

Consider the family of integrable Hamiltonian systems with Hamiltonian H{p, q\A) 
depending on the parameters A = (ai, . . . , cat). One can think of the classical Neu- 
mann problem of the motion on the unit sphere with the Hamiltonian 

H =^p^ + i(aigi H ^ aNq%) 

as an example. Let Fi{p, q\A), . . . , F^lp, q\A) be the corresponding N independent 
integrals in involution. 

Now suppose that we make the parameters time-dependent: A = A{t) with A{t) 
periodic with period T. In general this of course destroys integrability of the system. 
I was looking for the special cases, which I called integrable modulations, when the 
values of the integrals Fi{p{t),q{t)\A{t)), where {p{t),q{t)) is an arbitrary solution 
of the modulated system 

/ , N dp dH , , , , ss dq dH , , , , , - 

(1) ^^(^'''I^W)' 

are also T-periodic in time with the same period T. This implies that the time T- 
shift along the trajectories of the modulated system gives an integrable symplectic 
map with the same integrals i^i, . . . , Fn ■ 

Moser said that he would formulate the question as follows: when are the adia- 
batic invariants exact ? If we mean by the adiabatic invariants the action variables 
(see [T]), then this implies that the modulated dynamics preserves the Liouville 
tori foliation of the original integrable system, so we have simply non-uniform and 
non-straight motion on each of these tori. 

At that time I had an example of such integrable modulation for the rigid body 
system, which is the main subject of this section. But let me start with a simpler 
example of this phenomenon: the harmonic oscillatorQ 

If we consider the Hamiltonian of the harmonic oscillator in the form 

H=^{p^+LO^q^), 

then it is easy to see that it has no non-trivial integrable modulations at all. How- 
ever, if we consider the equivalent Hamiltonian system with 

(2) H^^cuip^ + q^), 
then any modulation lu = Lu{t) is integrable. 



I am grateful to Anatoly Neishtadt for pointing out this example to me. 
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Indeed, the modulated system 



dp dq 

clearly has integral 



- = -u.q, -=up 



(which is the action variable of the system), and by changing time t T{t) : ^ = 
Lu{t) we reduce it to the non- modulated form 

dp dq 
dr dr ^ 

The modulated system has the solutions 

p — RcosT{t), q{t) = RsmT{t), 

where the angle variable ip — T{t) depends on t in the way determined by io{t). Note 
that the Hamiltonian ([2]) (which is not an integral of motion anymore) changes 
periodically with time provided the modulation ^{t) is periodic. 

This example shows that in general for a prescribed dependence on the pa- 
rameters in the Hamiltonian we may not have integrable modulations at all. In 
particular, it seems that such modulations do not exist for the Neumann system or 
for the closely related Jacobi system of geodesies on ellipsoids [11 . That is why I 
found interesting that for the Euler rigid body system such modulation does exists 
even in the iV-dimensional version, which we are going to see now. 

After Arnold [l] the motion of an iV-dimensional rigid body is commonly iden- 
tified with the geodesic flow for the left-invariant metric on the Lie group SO{N) 

dX dX^ 
C = tr—J— . 

dt dt 

Here X G SO{N), J = is some positive symmetric matrix and denotes 
the transposition of matrix Y. The corresponding angular momentum AI e so{N)* 
satisfies the Euler equation 

(3) ^-[^^'"]' 

where il = X^^^ e so{N) is the angular velocity related to M by 

(4) M + .m. 

If we introduce the linear operator A : so{N) — > so{N)* by A{fl) = ilJ + Jft, then 
the Hamiltonian of the Euler equation can be written as 

H = ^trM'^A-^{M). 

The proof of integrability of the corresponding system was found by Manakov [S] . 
It is based on the important observation that the Euler equation can be rewritten 
in the Lax form 

(5) f^[L,P|, 

where 

(6) L^Af + AJ^ P^Q + XJ 
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and A is an arbitrary parameter (usually called spectral). Indeed, the relation ([5|) 
is equivalent to 

^ = [M, n], [j^,n] + [M, J] = 0, [J^ J] = 0. 

The third relation is trivial, the second one is equivalent to (|4]), so we end up with 
the original Euler equations. The Lax equation ^ implies that the coefficients of 
the characteristic polynomial of L 

P(A, fi) = det(L - ^/) = det(M + AJ^ - fil) 

are integrals of the systems. One can check that this gives enough integrals in 
involution to claim the Liouville integrability of the system and that the dynamics 
is linearisable on the Prym variety of the spectral curve 

Let me show now that Manakov's trick works also when J is not a constant but 
depends explicitly on t in such a way that 

with an arbitrary scalar function /(t), constant symmetric matrix Jq and the iden- 
tity matrix /. Indeed, the corresponding modulated Euler equation has the same 
form (l3|), but with t-dependent 

J=(Jo2 + /W/)i/2. 
The same arguments as above show that it can be rewritten as 

— = [M + xj\n + xj]. 

at 

For our particular J we have 

[M + AJ^ 17 + A J] = [M + XJ^ + f{t)XI, n + XJ] = [M + XJ^,n + AJ], 

which implies that the modulated Euler system can be written in the Lax form ^ 
with 

l^m + xjI, P^n + XJ. 

This means that the integrals of the non-perturbed system given by the coeffi- 
cients of the characteristic polynomial Pq{X, ^) = det(M -I- AJq — fil) are preserved 
by the modulated system, so 

J=(Jo2 + /(i)/)l/2 

is indeed an integrable modulation. The dynamics can be described explicitly as 
certain motion on the Prym variety of the corresponding spectral curve Po{X, fi) = 
0. 

Note that the Hamiltonian as well all other integrals of the initial system given 
by the coefficients of P{X, fi) = det{M + X,P — /i/) do depend on t: 

P{X,fi) = det{M + XJl + f{t)XI - fil) - Po(A,M - f{t)X) 

and their t-dependence will be periodic if the modulation f{t) was periodic. 

The corresponding T-shift gives an integrable discretisation of the Euler rigid 
body system, which is worthy of further investigation. It depends on the choice 
of f(t) and corresponds to certain shifts on the Prym varieties. In comparison 
with the discretisation from [T3] its generating function (Lagrangian) is probably a 
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very complicated function (cf. Bobenko-Lorbeer-Suris paper [3], where a different 
discretisation of the Eulcr top was proposed). 

3. INTEGRABILITY and ANALYSIS: IS THE CAT-MAP INTEGRABLE ? 

During one of my visits to Ziirich, when Sergei Tabachnikov was also a guest of 
FIM, we thought that it would be a good idea to interview Moser for the Russian 
journal for high-school students "Kvant" (American version "Quantum"), of which 
Sergei was one of the Editors. Moser agreed without hesitation. We gave him 
around a dozen questions, which we thought he would simply write the answers 
to. Instead he invited us to come to talk about this. I still can not forgive myself 
that we were not ready for this and not able to prepare anything for publication 
afterwards. Because to say that this was interesting is to say nothing. 

One of the questions was about mathematics and where he sees himself. Moser 
was very reluctant to divide mathematics into separate disciplines and considered 
himself simply as a mathematician. However he added that now he is not as active 
as before and left for himself a small garden, where he tries to follow all the main 
events. This garden was called "Analysis." 

Integrable systems is a crossroad area, naturally connected to geometry and 
algebra. Analysis is usually connected to it only through the KAM-theory. Jiirgen 
Moser is amongst the very few real analysts, who made crucial contribution to the 
modern theory of integrability. Many of us learnt this theory from his brilliant 
Bressanone's lectures jllj . 

I would like to talk now about one phenomenon in this theory which I have 
understood only recently. This is about the difference between smooth integrability 
and analytic integrability. One of the most famous Moser results is that KAM- 
theory works in smooth category as well as in analytic one, which was shown earlier 
by Kolmogorov and Arnold. Here is a simple example, inspired by Bolsinov and 
Taimanov which demonstrates how big is actually the gap between these two 
categories. 

Consider the following cat-map (or Anosov map) 



The classical Anosov result says that any hyperbolic torus automorphism is 
structurally stable in topology. This is the canonical example of a chaotic 
system with positive topological entropy, which in this case can be shown to be 
equal to log A, where A > 1 is the largest eigenvalue of A. 

I argue that there are several reasons to consider this system as integrable0 First 
of all, this is a linear map and therefore can be integrated explicitly by means of 
linear algebra. There is a problem with integrals though, because the dynamics is 

^To avoid complaints about abuse of terminology I would like to quote BirkhofT's classical book 
[2], whose 1966 edition was introduced by Moser: "When, however, one attempts to formulate a 
precise definition of integrability, many possibilities appear, each with a certain intrinsic theoretic 
interest." 



where 



T{? /T? is the two-dimensional torus and 
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ergodic on the torus. However I claim that the natural extension of cat-map to the 
cotangent bundle of the torus 

is integrable in the Liouville sense with smooth integrals. 

Indeed, in the eigen-coordinates u,v oi A the extended map has the form 

A : {u,v) (Am, A^^w), {pu,Pv) {X^^Pu, Xpv), 

so we have two integrals in involution: Fi = PuPv and 

=e"^sin(27r}^). 

log A"' 

Note that the second integral is smooth but not analytic. 

The original cat-map is sitting at the critical level Pu = Pv = 0- This is another 
lesson of this example (cf. Bolsinov-Taimanov [4J): Liouville- Arnold theorem does 
not tell us anything about what happens at the degenerate level of the integrals. 
In smooth category one might find a fully developed chaos there ! 

There is of course another, topological aspect of this example. Topological ob- 
structions to integrability have been discussed in the literature starting from the 
important work by Kozlov [8], but here I would like to make a different point. Let 
us look at the cat-map again. On the universal covering (which is the plane in the 
torus case) the dynamics is indeed integrable. Only after projection to the torus 
do we have ergodicity. A similar example is the geodesic flow on the surfaces of 
genus g > 1, which is chaotic on the surface but perfectly integrable on the hyper- 
bolic plane covering it. I would call such chaos topological. It has nothing to do 
with solvability. There are many examples of exactly solvable systems, which are 
ergodic. 

Let me discuss one of them (probably the oldest one) due to Claude Caspar 
Bachet (1581-1638), mathematician, poet, linguist, widely known by his transla- 
tion of Diophantus' "Arithmetica," where Fermat made his famous marginal notes. 
Looking for the rational solutions to the Diophantine equation 

y^-x^ = c, 

Bachet observed that if (x, y) is a solution, then so is 

This allows to construct infinitely many solutions to this equation. For example, 
for c = — 2 we have 

129 383 2340922881 113259286337292 
^' '^^ ^ ^100'" 1000^ ^ ^ 76602 ' 76603 > ^ •■■ 
The rational map B : CP^ CP^ 

— 8cx 

X B{x) 



4(0:3 + c) 

is called the Bachet map. Ceometrically this is the duplication map for the elliptic 
curve E given by y"^ ~ x^ — c. Indeed 

p(2^) = B{p{z)), 
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where p{z) is the correspondmg Weierstrass eUiptic function, so on E we have the 
map z ~f 2z. The Bachet map is actually the most chaotic among the rational 
maps in the sense that its Julia set coincides with the whole of CP^. In complex 
dynamics this important example of chaotic map is usually ascribed to Lattes, who 
worked at the beginning of the last century (see e.g. Milnor [10]). 

Again I would argue that this map has all the reasons to be considered as in- 
tegrable. One is exact solvability in terms of elliptic functions, which follows from 
the above arguments. Another is the existence of infinitely many commuting ra- 
tional maps (symmetries), which in the non-invertible maps replaces the integrals. 
Indeed, the rational maps _B„ defined by 

p{nz) = S„(p(z)) 

clearly commute with each other: 

BnBm = BmBn, B2 = B . 

Ritt showed that all such maps are related to some multiplication formula for elliptic 
functions. In many dimensions there are similar examples related to simple complex 
Lie algebras (see more details in [T7]). 

4. GeODESICS ON HYPERBOLOIDS: KnORRER MAP AND GEODESIC EQUIVALENCE 

AS "regularisations" 

Moser had a unique style as a lecturer. I attended one of his undergraduate 
courses on Dynamical Systems at ETH (initially to learn more German, but then 
just because I enjoyed it so much). Usually in the first hour he explained the result 
and the main ideas behind the proof using the simplest possible arguments and 
examples, then after the break full proofs with all the details were given. Not a 
single question was left unanswered and Moser, who usually kept time quite strictly, 
was always generous, making sure that that the matter was clarified completely. 
This combination of simplicity and clarity with mathematical honesty, not allowing 
him to hide the unwanted technicalities under the carpet, I found remarkable. 

The same concerns his writing. Moser told me that the large number (333) of 
derivatives in his result on invariant curves in KAM theory was required just for 
the sake of simplicity of the arguments and he could easily reduce it (I think he 
said to 18) at the cost of extra technicalities. The paper [12] is a nice example of 
Moser's remarkable style, which is relevant to the problem I would like to discuss. 

In the first part of this paper Moser considered Kepler's problem with Hamil- 
tonian 

By changing the time variable t to 




and using stereographic projection, he showed that the energy surface H = c with 
negative c can be compactified and that the regularised flow is equivalent simply 
to the geodesic flow on a unit sphere (three-dimensional in the usual case, when 
q G R^). We note that because of the collision orbits (corresponding to [gj ^ 0) the 
regularisation is really needed. The change of time allows to continue the motion 
after the collision in a smooth way. Moser mentioned that such regularisation was 
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essentially known to Levi-Civita and Sundman, but it looks like that this topological 
picture was clearly explained first in his paper. 

A simple remark I would like to make is related to the geodesic problem on 
hyperboloids. All the algebraic manipulations are the same as in the well-studied 
Jacobi geodesic problem on ellipsoids. However the behaviour of geodesies in these 
two problems are quite different: in the Jacobi case the dynamics is quasi-periodic 
while on hyperboloids the geodesies in general go to infinity. I would like to explain 
here that this difference disappears if one makes a "compactification-regularisation" 
of the geodesic problem on hyperboloid using Knorrer's transformation [7|. 

Let me recall first Knorrer's result (which I also learnt from Moser [H]). Let 

(7) {A-^x,x)^l 

be the equation of a quadric Q in the [n + l)-dimensional Euclidean vector space 
]R"+^, which for the beginning will be assumed to be an ellipsoid. Geodesies on Q 
in the natural (length) parameter s are described by the equation 

x" — XBx, 

where ' denotes the derivative with respect to s, i? = , A is the Lagrange 
multiplier determined by the constraint ([7|: 

{Bx, x') = 0, {Ex, x") + {Bx, x') = 0, \{Bx, Bx) -f {Bx\ x) = 0, 

Following Knorrer [7] , consider the change of parametrisation 
(9) s^r^ja[s)ds. 

where — —A. Knorrer observed the remarkable fact that the normal vector to 
the geodesic 

(10) 

\Bx\ 

in the new time r satisfies the classical Neumann system: 

(11) ^^^Bq + ^Iq, qeW^+\\q\ = l. 

In other words, the Gauss map (|10[) in combination with re-parametrisation Q 
transforms the geodesies into the trajectories of the Neumann system on the unit 
sphere in the harmonic field with the potential 

Uiq) = ^{Bq,q). 

The corresponding orbits have one of the integrals fixed: in Moser's notations [13] 

(12) *o(^,g)=0, 
where 

^o{u,v) = (1 - {Au,u))A{v,v) + (Au,vf,A = B^^. 



■^On the quantum level this corresponds to the 50(4) symmetry of the Kepler problem with 
negative energy discovered by V. Fock in 1935. 
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Thus, Knorrer transformation gives a bijection between the unit (co)tangent bundle 
Al to ehipsoid Q and the subvariety M of the (co)tangent bundle to the unit sphere 
given by relation (|12p . mapping the geodesies to Neumann orbits. 

Let us look now at the case when the quadric Q is a hyperboloid. In that case 
generically the geodesies will go to infinity and the corresponding A tends to as 
s — > oo. To see the rate of decay we note that the geodesic problem has the following 
Joachimsthal integral 

F = \Bx\^{Bx' ,x'). 

Indeed, 

F' = 2{Bx', Bx){Bx', x') + 2\Bx\'^{Bx', x") 
= 2{Bx', Bx){Bx\ x') + 2\\Bx\^{Bx', Bx) = 2{Bx', Bx){{Bx\ x') + X\Bx\'^) = 0. 
Using this integral we can write 

{Bx',x') _ F 



A = -- 



iBxP \Bx 



4 ■ 



Because at infinity \Bx\ sa cs we see that A w Cs"** and thus a = \/^A decays as 
(one should be careful with the sign of A in the hyperboloid case, but this can 
be done in all the cases). This means that the integral r = J a{s)ds is convergent 
at infinity, so it takes finite time in the new time variable r to reach infinity. The 
Neumann system tells us what happens afterwards: the geodesic will return from 
the "same infinity" (in general, along a different path), then go to a different one 
and so on. Note that the Neumann system does not change when B ^ B — zl for 
any z, so it does not feel the change from ellipsoid to hyperboloid, which affects 
only the value of integral ^'o (cf. Moser's comments in section 3.5 in |13)). 

The Knorrer transformation makes the necessary compactification and defines a 
smooth dynamics there. To see what happens geometrically, it is convenient to use 
the dual quadric Q* given by 

{Ay,y) = {B-'y,y) ^ 1, y = Bx e M"+i. 

Then the Gauss map becomes simply the central projection of the hyperboloid Q* 
to a unit sphere: 

y 
\v\ 

The image is an open set with the boundary B given by the intersection of the sphere 
with the asymptotic cone (Ag, q) — 0. The points at the boundary correspond to 
the "infinities" of the hyperboloid. Note that when {Aq,q) = the relation IT^ 
reduces to 

which means that the velocity vector must be tangent to the boundary. 
We can see this already in the simplest example of hyperbola 

b^xl + bixl — 1 

with 5o > 0, hi < 0. In this case the regularised fiow will be simply periodic back 
and forth motion on one branch of the hyperbola. The unit tangent bundle M 
in this case is topologically equivalent to the disjoint union of two open intervals, 
while on the Neumann side we have a circle, which is a compactification of M. with 
two extra points, corresponding to two "infinities" of the hyperbola. 
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The two-dimensional case is more interesting. In particular, in the case of a one- 
sheeted hyperboloid we have two families of straight lines, which are the simplest 
geodesies. For them (and only for them) we have {Bx' , x') = and correspondingly 
A = 0, which means that the Knorrer transformation formally can not be applied. 
The central projection maps them into large semicircles, corresponding to the large 
energy limit of the Neumann system. This shows that we may have something 
extra on the geodesic side as well. 

There is another, more geometrical way to compactify the geodesic system on 
hyperboloid, which I have realised after discussion of this problem with Alexey 
Bolsinov. It is based on the notion of geodesic equivalence of the metrics. Two met- 
rics ds\ and ds\ (Riemannian or pseudo-Riemannian ) on a manifold M" are called 
geodesically equivalent if they have the same geodesies considered as unparametrised 
curves. A trivial example is given by the proportional metrics ds\ — cds\, but there 
are other examples as well. The problem of description of projectively equivalent 
metrics has a long history, going back to the work of Beltrami and Dini in XlX-th 
century (see the recent review [16j . where the role of integrability in this problem 
is emphasized). 

We need the following example of two geodesically equivalent metrics, which, 
in spite of its classical look, seems to be discovered only recently independently 
by Tabachnikov [13] and Matveev-Topalov (TB]. Let ds\ be the restriction of the 
Euclidean metric 



ds^ = dxl — {dx, dx) 



i=l 



in R"+^ to the quadric Q given by ^ (which is a standard metric on Q) and 
be the restriction of the metric 

fl3) dr' - Mxf +Ma;i + --- + &„+irf4+i _ _J_(Bdx dx) 

to Q. Then the claim is that these two metrics are geodesically equivalent (see 
[15l[T6]). This result was originally formulated for the ellipsoids, but it obviously 
holds for the hyperboloids as well although the second metric may become non- 
positive (cf. Theorem 4.13 in Khesin- Tabachnikov [6]). 

The observation, which may possibly be new, is that in the hyperboloid case the 
second metric ds^ can be naturally extended to a regular metric on the projective 
closure Q C MP"+^. The easiest way to see this is by direct calculation. In the 
projective coordinates : zi : ■ ■ ■ : Zn+i the quadric Q is given by the equation 

bizl H h6„+iz^+i = zg, 

where the original quadric Q is the affine part in the chart, where zq 7^ and 
Xk = Zk/ zq. In a different chart, where say zi ^ 0, we can use the affine coordinates 
Vk = Zk/ zi, fc = 0, 2, . . . , n -f 1 such that 

_ J_ _ y2 _ Vn+l 

X\ — , X2 — 7 ■ ■ ■ 7 '^n+l — 7 

yo yo yo 

in which the equation of Q becomes 

61 -I- 621/2 ^ bn+iyl+i = yo- 
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One can easily check that the second metric cis| in these coordinates coincides with 
the restriction to Q of the metric 

(note that outside Q the metrics dr^ and dr^ are different). A remarkable fact 
is that this metric is regular when j/o 0, which means that ds^ indeed can be 
regularly extended to the whole Q. 

When n = 2 we have one-sheeted and two-sheeted hyperboloids when 6i < < 
&2 < &3 and 6i < 62 < < 63 respectively. In the first case Q — topologically 
is two-dimensional torus and the metric ds^ has the signature (1, 1), in the second 
case Q = S"^ is a topological sphere with the Riemannian metric {—ds^). 

The relationship of the geodesic equivalence with Knorrer transformation and 
the general geometrical picture is worthy of further study. We only mention here 
that the comparison of the formulas ([8]) , ^ with ([T3|) shows that the corresponding 
reparametrisations are the same. It may give also an alternative explanation of the 
special role of metric (|13|) on the quadrics. 

Another interesting question is to look at the properties of the spectral problem 

-V'" + E\{s)i) = 0, 

where A is the corresponding Lagrange multiplier ^ (see OHH], where the ellipsoid 
case was considered). 
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